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Abstract 



We present a new measurement of the branching ratio of the Cabibbo suppressed 
decay D° — > 7r~fj, + u relative to the Cabibbo favored decay D° — > K~ ji + v and an 



improved measurement of the ratio 



/+(Q) 
1W) 



Our results are 0.074 ± 0.008 ± 0.007 



for the branching ratio and 0.85 ± 0.04 ± 0.04 ± 0.01 for the form factor ratio, 
respectively. 



1 Introduction 



Semileptonic decays provide the advantage of having factorizable weak cur- 
rents in the Hamiltonian which allows for a clean theoretical description. 
The hadronic current can be described in terms of two form factors, f+(q 2 ) 
and f~{q 2 ) which are functions only of the lepton-neutrino invariant mass 
squared, q 2 . Assuming a pole dominance parameterization of the form fac- 
tors, we present a parametric analysis of the pseudoscalar semileptonic decays 
D° — > 7r~ ji + h> and D° — > K~fi + v from the FOCUS experiment. 

This paper concentrates on the relative branching ratio and the form factor 
ratio of the Cabibbo suppressed decay relative to the Cabibbo favored mode. 
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Since the efficiency tends to have a non-negligible q 2 dependence (see Figure 3), 
we allow the pole masses and the ratio /_(0)// + (0) to vary freely in the fit. 
The results and description of the detailed analysis of the pole masses and 
/_(0)//+(0) are included in another paper [1] along with a non-parametric 
analysis of the high statistics decay D° — > K~fi + v . 

We report a new measurement for the branching ratio F(Tr~ fj, + u) /F(K~ fj, + u) 
in agreement with recent results from the CLEO collaboration [2,3]. These 
results indicate a lower value for this branching ratio than the one reported 
in the PDG [4]. We also report a new measurement of the form factor ratio 

/ ir (o) 

jtj^ with greatly improved errors with respect to existing measurements 

and compare it to recent theoretical predictions from an unquenched Lattice 
QCD calculation [5,6] . 



2 Data Selection 



This analysis is based on data collected by the FOCUS experiment during 
the 1996-97 fixed target run at Fermilab. FOCUS is a photoproduction ex- 
periment which collected a large sample of charm decays produced in the 
interactions of a photon beam [7] with an average energy of ~ 180 GeV on 
a BeO segmented target. The FOCUS spectrometer [8,9,10,11] is equipped 
with a 16 plane silicon strip vertex detector; 4 planes are interleaved with the 
targets and 12 planes are located downstream of the target area. Momentum 
analysis is accomplished by two magnets with opposite polarities and 5 mul- 
tiwire proportional chambers. Three multi-cell threshold Cerenkov counters 
provide charged particle identification. A muon counter located at the end of 
the spectrometer is responsible for muon identification. 

We reconstruct the semileptonic decays D° — > 7r~ fi + u and D° — > K~fi + v re- 
quiring a D* + -tag where the D* + is reconstructed in the D°tt + final state. 1 
Whenever possible we apply identical selection criteria to both decay modes 
to reduce systematic effects. As the decay D° — > n~fi + is has more background 
and less statistics, the selection cuts have been optimized for this mode. The 
signal and normalization samples are selected requiring two opposite charged 
tracks form a good vertex with a confidence level greater than 1%. One of 
the two tracks from the D° decay vertex must be identified as a muon from 
the inner muon detector with a confidence level greater than 1% and must 
have momentum greater than 10 GeV/c. To suppress pion and kaon in-flight 
decays, this track is required to have a consistent momentum when measured 
in the first and second magnets. The other track must satisfy a Cerenkov re- 



Throughout this paper charge conjugate modes are implied. 
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quirement based on the value of the negative log-likelihood W for a given 
hypothesis: in the D° — > 7r - /i + z/ mode, the pion must be favored with respect 
to the kaon hypothesis by at least 3 units of likelihoods (W(K) — W(ir) > 3); 
in the case of K~/j, + u the kaon must be favored over the pion hypothesis by 3 
units of likelihoods (W(n) — W(K) > 3). To reduce non-charm background, 
the candidate hadron must have a momentum greater than 14 GeV/c. The 
primary vertex is found after excluding the candidate tracks from the D° 
decay vertex; the remaining tracks are used to form candidate vertices. Of 
these vertices we choose the one with highest multiplicity and we break am- 
biguities by picking the most upstream vertex as the primary vertex. This 
vertex is required to be isolated from other tracks in the silicon strip ver- 
tex detector by requiring that the confidence level of any another track not 
used in the determination of the primary or decay vertex be less than 1%. 
For each hadron-lepton combination that satisfies the above requirements, 
another track coming from the primary vertex must be found as the candi- 
date "soft" pion from the D* + — > D°tt+. The 7r+ candidate must not have the 
pion hypothesis strongly disfavored over all other particle hypotheses from the 
Cerenkov system (min{W^(e), W(K), W(p)} — W(ir) > —6). It must also have 
a momentum greater than 2.5 GeV/c. To suppress backgrounds from decays 
where a final state particle is lost (usually ir°), such as K~n + -K°, K~it + it it , 
p~fi + u and Kirfi + i> 2 , we place a lower cut on the hadron-lepton invariant 
mass (visible mass) of 1.0 GeV/c 2 . Contamination from D° — > K~ir + is elim- 
inated by requiring the visible mass to be less than 1.7 GeV/c 2 . 

Since the neutrino is not reconstructed, the resultant smearing effects on the 
resolution play an important factor in this analysis. Rather than using the 
standard neutrino closure resulting in a two-fold ambiguity on the D° momen- 
tum, we take advantage of the -D* + -tag by boosting the final state particles in 
the hadron-lepton center of mass frame. By constraining the K~ ji + u ji + u ) 
mass to the D° mass and the K~/j, + u 7r+ (tt~ij, + u 7r+ ) mass to the D* + mass, 
we are able to determine the angle between the neutrino 3 and the 7r+ direc- 
tion. We then sample the azimuthal angle and choose the one that gives the 
direction of the D° most consistent with pointing to the primary vertex. 



3 Analysis 

The fit to the data is designed to constrain the background in the ir~ fi + u sam- 
ple and to supply information about the pole mass and form factors. To accom- 

2 With the notation Rtt^u we refer to the sum of the decays to K°7r~ p + v and 
K~ir° ' n + v from D° or to the sum of the decays to K~ir + n + v and K°7r°// + i/ from 
D+. 

3 The neutrino and the D° directions are the same in this reference frame. 
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Fig. 1. D* + — D° mass difference distributions for D° — > TT^fi + u (left) and 
D — > K~ n + v (right). The amount of non-peaking background is found from a 
fit to q 2 and D* + — D° mass difference distributions. The vertical line indicates 
where the cut is placed. 

plish these goals we perform fits on two-dimensional distributions where the 
free parameters are the signal and background yields. All the fits are binned 
maximum likelihood fit where the likelihood is defined as: 



where (n» 3 -) is the number of expected (observed) events in the bin ij. 
First, a fit of q 2 and D* + — D° mass difference is performed to establish the 
amount of non-peaking background (Fig. 1). 4 We next place a mass cut on the 
D* + — D° mass difference of less than 0.154 GeV/c 2 to reduce the background 
and to obtain more reliable results for parameters such as pole masses and 
form factors. A fit is then made to the two-dimensional distribution q 2 vs. 
cos 6e (where cos 9e is defined as the cosine of the angle between the neutrino 
direction and the D° direction in the rest frame of the lepton- neutrino system). 
The fit is first performed on the K~fi + u sample and the results from this fit are 
used to set the background from K~ /i + z/ and K-n^v in the 7r~/i + z/ sample. 

4 We define "peaking background" in the ir~ sample as the sum of the back- 
ground contributions from D° — > K~ [i + v , K-k^l + v and p~ [i + v , while in the K~ fi + v 
sample the peaking background is given only by K~ifi\i r v. 
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In the fit to the K~ fi + u distribution we make use of the recent vector to pseu- 
doscalar branching ratio measurement T(D + K^p + v)/T{D + — > K°ii + v) = 
0.63 ± 0.05 [12] in the form of a penalty term added to the log-likelihood as 
shown in Eq. 2: 



0.63 



. Y„- + e(K-Kix+v) 

/\,, = 2 lo K £ Klll , + v " + » ( (Q ^ (2) 



where we assume isospin invariance to relate D + and D° decays. The likelihood 
C is constructed using the expected number of events is each ij bin of the two- 
dimensional distribution given by: 



where in Eqs. 2 and 3 the fit parameters Y a are the fitted yields, S a are 
the normalized shapes obtained from Monte Carlo and e the reconstruction 
efficiency We define the cc component as the background obtained from a 
high statistics charm-charmbar Monte Carlo sample after removing the modes 
handled specifically in Eq. 3 (and 5). 

In a similar way we fit the 7r~fi + v distribution. We use the branching ratio 
T(D° -> p-/i + is)/T(D° -> fTvr/i+z/) = 0.086 ± 0.010 5 to constrain the back- 
ground in the fit as shown in Eq. 4: 



Y p -v+u e(Kn»+v) _ Q QgfiV 

TP - _9 W C 4- V r ^+. <P " + ") )_ ( ^ 

■ ' -//// - - ■",-> —:■!■■" (0 010) 2 ' 



The expected number of events in each two-dimensional bin used to construct 
the likelihood is: 



5 We estimated the branching ratio of D° — > p~ p + v relative to D° — > Ktt[i + v using 
the weighted average of a recent result from the CLEO-c collaboration [3], the PDG 
values [4] and a preliminary result from FOCUS [13] where we correct for isospin 
when necessary. 
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"T r A"-ir°/x + " 



e((X- -> 7r-)7rV + ^) 



:' l J 



(K-^nr-)ir n+v~' r 



(5) 



where Y^- +ll and Y^ w0 +l/ in Eq. 5 are fixed to the results obtained from the 
fit to the K~/j, + i> data (Eq. 3). The symbol (X — > Y) means that a hadron X 
is misidentified as Y. 

To measure pole masses and the form factor ratio rj = /5(0)//^(0) we apply 
an event-by-event weighting procedure [14]. This is achieved by re-weighting 
each Monte Carlo event according to the ratio of the probability that the event 
was generated with a pole mass M' ole and a form factor ratio rf relative to the 
probability that the event was generated with the default values M D * (M D * 
for nfiu ) and r]°. 6 The relative efficiencies of the decays D° — > ir~fi + v and 
D° — > K~/jl + u are defined as the ratio of the reconstructed and generated 
Monte Carlo events. At each fit iteration these efficiencies change as a function 
of the pole masses and rj values. 

The weight Wi for an event with q 2 = q 2 is given by the equation: 



= I(M polc , V ';q 2 ) N(M D ., V °) 
1 I{M m ^-q 2 ) N(M^rf) U 



where the intensity is: 



J(M pole , m q 2 ) oc fl(M polc ; q 2 ) g(rj) (7) 
and the normalization is determined by: 



Ngen 

N(M polc , V ) = £ f 2 + (M pole ;q 2 ) g( V ). 

i=l 



6 The default values for the parameter r] are: if = —0.724 for D° — > i^T /i + z/ and 
r? = -0.856 for D° -» 7i-/i+i/ • 



7 




(GeV/c 2 ) 2 



-0.5 n 
COS0 



(GeV/c 2 ) 2 



Fig. 2. Fit projections for -k~/j, + i/ and D° — > K~/j 1 + u. The fit is performed on 
a two-dimensional distribution of g 2 and cos 0g. In the D° — »■ tt~^ + u, the peak- 
ing background contribution is defined as the sum of the contributions from 
D° — > K~n + v, p~ [i + v and K7t/j, + v. 

The form factor / + (M pole ; g 2 ) is assumed to have the following q 2 dependence: 



/+(M pole ;g 2 ) 



/ + (0) 
1 - 



(9) 



M 



pole 



and g{rj) can be written in terms of three kinematic coefficients A, B and C: 
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g{rj) = A + B r} + C r} 2 . 



(10) 



From the fit to the 7r _ /i + z/ (K~fi + v) distributions (Fig. 2) we find 288 ± 29 
7T~fi + v (6574 ± 92 K~fi + v) events. Correcting for the relative Monte Carlo 
efficiency we find the branching ratio for the Cabibbo suppressed decay D° — > 



7 The kinematic dependence is shown in detail for kaon semileptonic decays in 
reference [4] on page 618. 
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Fig. 3. Reconstruction efficiency as a function of the q 2 for -k~{i + v (top) and K~ fi + v 
(bottom). 

7r~ relative to the Cabibbo allowed decay D° — > K~fi + v to be: 



r(D°- 



■7T /i Z/J 



r(D°^ir-/i+^) 



0.074 ±0.008 (stat.) 



1.93±8;8| for the 



From the same fits we find M w = LQlTo^ and M K 
TT~fi + v and the K~ pole masses respectively. We also measure the ra- 
tio /5(0)//£(0) = — l-7±i;|. A detailed description of the pole mass results 
has been included in Ref. [1]. 

Using the yields from the fit it is possible to obtain the ratio of the form 
factors /+(0)//^(0). In order to do this we compute a numerical integration 
of the differential decay rate modulated by the reconstruction efficiency as a 
function of the q 2 [15]. This efficiency is found by sampling the q 2 Monte Carlo 
distribution and dividing the reconstructed events by the generated events in 
each bin. The resultant distribution is then fit to a third degree polynomial 
(Fig. 3) which is used in the computation of the integral. We quote the result: 



■fflO) 

mo) 



0.037 ±0.004 (stat. 
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Applying the unitarity constraints on the CKM matrix elements [4] we use 
the value |^| 2 = 0.051 ± 0.001 in Eq. 12 and measure the ratio /^(0)//^(0) 
to be: 



= 0.85 ±0.04 (stat.). (13) 



4 Systematic Studies 

Several studies have been performed to search for possible systematic uncer- 
tainties. The fitting procedure was tested on a Monte Carlo set whose size is 
roughly 20 times the FOCUS data set and we verified that the fit returned 
the input values used in our simulation. 

We checked for possible biases as well as the accuracy of our statistical error by 
performing a fit on fluctuated data distributions multiple times and comparing 
the mean and width of the distribution of the fit results to our measurement. 
We found that we have to add a 0.005 contribution to the systematic error 
to compensate for K~ n + u and Kn^u contributions that were not allowed 
to float in the 7i'fi + u fit. We also performed an analogous study using the 
fit function as the parent distribution to establish how well our fit function 
described the data. We compared the likelihood obtained from our measure- 
ment to the distribution of the likelihoods from the fluctuated fit function. 
We found good agreement indicating that our fit function well represents the 
data. 

We investigated the stability of our results by changing a variety of selection 
criteria: the significance of separation between the primary and secondary ver- 
tex, muon identification, track momenta, visible mass cut, and Cerenkov iden- 
tification. We found no significant change in our results and assign a systematic 
uncertainty of 0.003 on the branching ratio due to cut variations. This number 
is found by computing the variance of this set of results. 

We further investigated fit variations by using a different approach in which 
we fit the q 2 and D* + — D° mass difference. Rather than fitting the K~ n + u 
distribution first, this fit was performed simultaneously on the ii~ fi + u and 
K~fi + u samples. The results are nearly identical to the results obtained from 
the fit to q 2 and cos 9e- Other fit variations include changing the bin size. By 
computing the variance of these a priori likely results, we assigned a systematic 
uncertainty of 0.004 from fit variations. 

Since the Monte Carlo is used to determine the amount of K^fi + u background 
in the n~ ji + u sample, we are sensitive to the simulated misidentification rate. 



mo) 
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We used the high statistics modes D° — > K~n + and D° — > fT~7r + 7r + where no 
Cerenkov requirement was applied to measure the K ^ n misidentification 
rate and we used the statistical error on the combined sample (after applying 
the same Cerenkov requirement used to select the D° — > n~fi + is events) to as- 
sign a systematic uncertainty. We varied the misidentification rate so obtained 
by ±lcr, and we find a contribution of 0.002 to the total systematic error. 

The contributions to the systematic error on the branching ratio and the 
corresponding contributions to the error on the form factor ratio /£(0)//f (0) 
are listed in Table 1. 





BR 




Fluctuated data distribution 


0.005 


0.029 


Cut variations 


0.003 


0.017 


Fit variations 


0.004 


0.023 


Cerenkov misidentification 


0.002 


0.011 


Total 


0.007 


0.042 



Table 1 

Sources of systematic errors and relative uncertainties. The contributions to the 
error on the ratio f+(0)/ f+ (0) are found by propagating the corresponding errors 
on the branching ratio. 

In the measurement of the form factor ratio f+(0)/f+(0) we also added vari- 
ations on the fit to the reconstruction efficiency as a function of the q 2 used 
in the numerical integration. We varied the bin size and fit functions. We find 
a contribution to the systematic error of 0.010 which is added in quadrature 
to the errors propagated from the branching ratio measurement. 



5 Summary and Conclusions 



We quote the final results as: 
r(£>°^7r>+z/) 



0.074 ± 0.008 (stat.) ± 0.007 (sys.) 



(14) 



and 



cd 



Vc 
Vr, 



0.037 ± 0.004 (stat.) ± 0.004 (sys.). 



(15) 
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Using \y i \ 2 = 0.051 ±0.001 from unitarity constraints, we find the form factor 
ratio to be: 



mo) 



0.85 ± 0.04 (stat.) ± 0.04 (sys.) ± 0.01 (CKM) 



(16) 



where the last error (CKM) corresponds to the uncertainty on the ratio |lcd/Ks| 
We compare our results to the measurement reported by the CLEO collab- 
oration in Ref. [2] where they report the branching ratio of D° — > 7r~e + z/ 
relative to D° -> K~e + v to be 0.082 ± 0.006 ± 0.005 and the form factor ratio 

jw^ = 0.86 ± 0.07^q;o4 ± 0.01. We also compare our branching ratio result 

;o the recent measurement from absolute branching ratios for D° — > 7r~e + z/ 
and D° — > K~e + v from CLEO-c [3] where they report a relative branching 

ratio of 0.070 ± 0.007 ± 0.003. Our results are consistent with both of these 

mo) 



n (o) 



new measurements. Further, we report an improved measurement of 

in good agreement with SU(3) breaking expected in recent lattice QCD cal- 
culations where they quote a form factor ratio value of 0.85 ± 0.05 [5] and 
0.86 ±0.05 ±0.11 [6]. 
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